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In  this  study,  cocvllcaied  mchine  plant  Interactions  were 
analyzed  so  that  sore  efllclent  harvesters  could  be  designed.  Plants 


only  where  a ehealcal  or  biological  method  is  prohibited  or 


CHAPTER  1 
INTRODUCTlOtl 


emergent  hydrophytes  or,  in  neeh&nlcal  terms,  deFornible  bodies  that 
are  positively  buoyant.  These  plants  grow  together  Co  Torn  a nit  or 


{Schardt,  1986).  Florida's  Surfaeo  Voter  leprovement  and  Honogemeri 
Act  (1987)  Includes  aquatic  plant  mnogononl  strategies  to  restore 


plants  and  anioals  (Rosiey,  198T).  'These  political  issues  and 


1967).  'Typically.  200 


tons  of  plants  are  found  In  one  hectare:  consequently, 

createenc  costs  approxioBtely  SlOO/ha.  The  I 
Enslneere  considers  a systn  which  harvests  and  disposes  of  60  to  lOO 

the  water  hyacinth  (Culpepper  end  Decell,  1976).  The  Corps  also  fount 

conditions,  and  that  the  principal  problen  Is  the  transportation  of 

towing  a eat  of  plants.  Often,  however,  when  e sat  is  towed,  the 
plants  lend  to  "Jam",  Bogart  (1949)  described  the  ntechanlsins  of  Jtus 
fortsatlon  as  a cooibinatlon  of  accretion  of  plants  and  a rolling  under 

consequence  of  n Jam,  towing  force  increases,  ond  occoslonally  a eat 

equipment  that  Is  based  on  hydromechanical  characteristics,  (e,g,  draj 

physical  characterization  of  aquatic  plants  has  been  inadequate  to 
predict  the  performance  of  large  towing  or  removal  systems.  The 
purpose  of  this  research  was  to  determine  hydromechonlcal  properties 


ajAFTER  II 
OBJECTIVES 


The  principle  objectives  ol  this  research  are  listed  beloe.ond 
each  objective  will  be  Che  focus  of  a subsequent  ct^pcer. 

-Ihuuititatlvely  describe  chose  physics]  characteristics  of  water 

-DeCernlne  Che  hydrodynaiilc  and  elastic  contact  properties  of 
water  hyacinths,  and  find  the  equations  chat  describe  the 

equations,  a imthewBClcal  model  that  predicts  how  a particular 
a^gre^clon  of  water  hyacinths  reacts  to  an  applied  velocity. 

-Validate  the  model  and  demonstrate  practical  applications  for 


CHAFFER  III 

WYSICAL  CHABACrausriCS  OF  fATER  HYACINTH 
ACCRECATIOnS 


Several  physical  characterlscies  of  individual  water  hyacinth 
plants  and  of  different  a^regacions  are  described  la  this  clmptcr: 
standing  density,  plant  mss  and  structure,  sac  buoyancy  and 
connectivity.  The  latter  is  a new  property.  Generally,  these 
characteristics  vary  according  to  habitat,  season  and  plant  age,  and 
affect  harvest  operations. 


li  description  of  m 


provides  expected  ranges  oF  plant  mss;  leaf,  root,  and  rhiacsse  length 
and  diameters,  and  the  internal  cellular  structure  of  the  plants. 

This  Infonsatlon  is  boalc  for  the  detemlnatlon  of  inertia  {e.g,  plant 
mass)  and  for  the  Identification  of  the  plant  parts  that  physically 


projecting  into  c 


principal  plant  parts 
iscQus  forces  generated  when 
to  predict  rhizome  lengths 


difficult 


they  ere  usually  heavily  covered 


Standing  blonaas  density,  In  combination  w: 
and  vidch,  determines  the  capacity  of  a harvester  and  the  power 
requlreoiencs  for  conveying  or  cranaportation  systems.  Also,  plant 
mass  affects  Inertial  forces  eaeouniered  in  a towing  operation.  On 
the  other  hand,  areal  population  density  is  useful  in  flow  models 
because  total  projected  area  In  flow  depends  on  the  number  of  plants. 
The  measurement  of  areal  population  density  Is  time  consuming  because 
It  requires  handling  every  plant;  yet,  an  empirical  method  to  predict 
population  density  from  other  physical  description  (e.g,  plant  mss  or 
length)  does  not  exist, 

Hat  buoyancy,  the  force  required  to  submerge  a given  area  of 

under  such  devices.  Water  hyacinth  plants  can  be  connected  through 
leaf  entanglements  and  stolons.  Hats  have  been  characterized  by 
either  standing  population  or  biosass  density  and.  to  a limited 
extent,  by  buoyancy,  but  the  literature  reports  no  property  that  would 
quantify  the  "connecclveness"  of  a mat;  therefore,  a test  described  in 
this  chapter  was  developed  to  measure  mat  connectivity-  Information 
on  connectivity  is  needed  when  a mot  is  to  be  separated  into  smaller 
units  for  towing  or  lifting  from  the  water-  Also,  information  on 
connectivity  changes  could  be  useful  In  biological  studies.  Rapid 
decreases  In  connectivity  in  a mat  could  Indicate  that  the  plants  have 
been  subjected  to  physiological  stresses. 

-Scale  frOTi  available  literature  the  structure  and  composition  of 


hyacinths. 


-Determine  the  distribution  of  mass  in  a vatei 
aggregation  and  hoe  that  distribution  could  vary  fi 

-Develop  a relation  between  standing  density  c 


hyacinth 

in  ana  population 


-Develop  an  empirical  method  to  calculate  rhizome  length. 
-Determine  how  mat  buoyancy  varies  for  different  rats,  and  find 
relation  with  the  use  of  plant  size  characteristics  to  predict  it. 

Review  of  Literature. 

ICorpho-anatonleal  Description  of  the  Woter  Hvoclnlh. 

The  nature  water  hyacinth  plant  is  a free-floating  (usually), 
emergent  hydrophyte  chat  consists  of  roots,  rhizomes,  stolons,  leaves 
Inflorescence  and  fruit  clusters  (Figure  3.0).  In  Uniting  nutrient 
conditions,  shoot  to  root  length  ratio  In  the  range  of  1-2  is  coBsnon, 
and  53X  of  the  total  plant  dry  matter  Is  In  the  leaves  and  petiole. 

In  nutrient-rich  conditions,  shoot  to  root  ratios  are  in  the  range  of 
d-6,  and  up  to  85*  of  the  dry  matter  Is  in  the  foliage  (Reddy,  1988). 
MediUfs-slzed  sacer  hyacinths  (500-2000  g/n^  db)  are  harvested  from 
sewage  treatment  plants  because  plant  densities  are  optimum  for 
achieving  maximum  biomass  yields  (Reddy.  1988). 

The  fibrous  root  (as  opposed  to  a tap  root)  has  a dense  oluster 

dispersed,  and  they  are  arranged  in  groups  alor«  the  rain  root  axes 
(Mallory  et  al..  1970).  The  roots  of  plants  from  different  habitats 
and  seasons  vary  little  In  diameter  but  can  vary  greatly  In  length. 
The  diameter  varies  from  O.d  Co  1.6  nil.  The  length  varies  from  4 to 


15  cm  In  siruU  pleuiCs,  10  to  36  cn  In  mediun  plants,  and 
In  largo  plants  [Ponfound  and  Earle,  1948  and  Zajnora  at  a 


The  cortical  regions  consist  of 


A leaf  oonstscs  of  o petiole,  a float  (In  snail  plants  or  those 
living  Id  very  open  areas)  a blade  and  possibly  an  Isthnus  (tenuous 
portion  between  the  float  and  the  blade).  Floats  are  located  In  the 
petiole.  They  are  round  or  pear-shaped  and  range  from  2 to  4 cm  In 
diameter.  They  are  formed  when  the  phyalologltal  availability  of 
water  to  the  plant  Is  great  (Roa,  1920):  hence.  In  very  dense  stands 
of  water  hyacinths,  floats  are  usually  absent.  Petiole  length 


with  Increasing  plant  standing  density  (Tucker,  1981). 


Iis/ii  (>b),  petiole  length  eas  22  cm:  at  20  kg  /a  (»b).  lei^th  was 
and  at  40  kg  /a  (wb),  length  was  50  cm-  Generally,  the  float 
typo  of  leaves  are  disposed  in  a nearly  horizonal  position  (IS  to  46* 
from  horizontal)  wheroes  the  equltant  leaves  (e-g.  leaves  whose  bases 
overlap  within  or  above  then)  approach  verticalicy  (75*  from  the 
horizontal)  (Penfound  and  Earle.  104B).  In  dense  stands,  water 
hyacinths'  leaves  grow  upright  rather  than  horizontally  over  the  eater 
surface  (Reddy,  1988). 

hyacinths  at  their  distal  ends.  Stolons  can  be  relatively  short  and 
vertical  in  dense  pats  lait  long  and  horizontal  In  open  conditions. 

5 to  8 cm  long  in  close  stands  (Penfound  and  Earle,  194S  and  Zamora  et 
al..  1382). 

Internal  structures  of  various  organs  of  the  aerial  part  of  a 
plant  are  similar  to  each  other  (Zamora  et  al..  1982),  and  they  can  be 
classified  into  three  groups:  parenchymatous  cells,  support  and 
conduction  cells,  and  air  spaces,  fiinety-flve  percent  of  the  bloeass 
of  a plant  is  parenchymatous  or  consists  of  unspoclallzcd  cells  that 
usually  have  well-defined  central  vacuoles,  and  they  are  used 
typically  for  storage  of  protein,  starch  grains,  oil,  etc.,  or  as 
chlorenchyna  cells  In  the  blade. 

Collateral  vascular  bundles,  a collection  of  phloem  and  xylen 
cells,  are  scattered  throughout  the  aerial  part  of  a plant.  Phloem 
functions  in  conduction  of  food,  and  xylem  functions  in  conduction  of 

and  trachelds,  which  are  types  of  xylem  providing  mechanical  support. 


blades,  roots  and  petioles,  but  rare  In  stolons  and  rhizoiies,  while 
tracbeids  can  be  located  throughout  an  entire  plant.  The  fibers  are 
1.5  to  3 nr  long  by  2 to  10  pn  in  dlojDcter,  whereas  trachelds  are 
elongated  and  helical  In  shape  and  can  be  stretched  far  before  they 


Air  spaces  allowing  water  hyacinths  to  float  are  ioeatod  In 
abundant  Intercellular  spaces,  air  chaoibers  and  parenchyma  cells.  Air 
chambers  or  lacunae  are  particularly  characteristic  to  perenchyim 

been  pulled  opart  at  the  corners  by  the  continuous  growth  of  ths 
epidermis.  The  air  chambers  formed  permits  diffusion  of  gases  and 
equalized  pressure.  The  cross  portions  of  the  air  chambers,  known  as 
dlaphragre,  prevent  flooding  because  those  structures  are  perforated 
with  minute  openings  specific  only  for  the  passage  of  gases.  The 


Basically,  the  water  hyacinth  consists  of  water.  Penfound  and 
Earlo  (1948)  and  Knlpllng  et  al.  (1971)  determined  that  93  to  95*  of 
cho  fresh  content  of  roots,  95. IX  of  rhizomes,  96. 7X  of  stolons,  93. 9X 
of  floats,  and  67  to  B9.3X  of  blades  Is  water. 

Depending  on  the  standing  density  of  plants,  50  to  65X  of  the  dry 
weight  of  a plant  Is  fiber  (cellulose,  hemicellulose,  and  lignin), 
around  40X  is  carbon,  and  1.3  to  3X  is  nitrogen,  (Oebusk  and  Dierber. 


1984). 


irle  (1948).  He  found  chat  i 


specific  eeight  of  a root  was  0.782:  of  a rl 


■hen  othcntlse  noted,  seven  plants  were 
coeiponent  description  per  nal  or 


Plant  Hess  Distribution  and  Rhlloise  Length. 


previous  section  except  that  nire  than  ten  plants  irere  sampled  (see 
results  for  exact  numbers)  and  rhizomes  were  measured.  Itats  from 
these  locations  sere  sampled  to  obtain  data  that  represents  different 
plant  sizes  and  conditions.  The  plants  from  rhs  two  lakes  and  from 
Zellwood  had  long  roots.  The  plants  from  Bivan's  Am  were  very  big. 
Swine  unit  cultures  were  harvested  periodically,  and  they  had  short 

The  classification  of  plant  size  within  zats  was  attempted  by 
this  method:  first,  sample  populations  of  plants  were  measured  to 

determine  if  the  mass  data  was  distributed  norimliy  (i.e.  skewness  was 
calculated).  Once  this  was  ascertained,  the  mass  values  were  attached 
Co  plant  size  classes  and  their  frequency  noted.  A histogram  for  each 
sample  population  of  water  hyacinths  was  plotted  using  the  classes. 
Skewness  was  calculated  from  each  frequency  table  to  check  if  the 
diecrlbuclon  of  the  size  classes  within  aach  ntts  was  norml. 
ConnetHvlty,  Buoyancy,  and  Areal  Densities. 

After  a physical  analysis  of  water  hyacinth  mats,  it  was  observed 
that  mats  could  differ  by  average  size  of  plants,  standing  density, 
eat  buoyancy  and  connectivity.  Except  for  connectivity,  there  were 
established  methods  to  measure  all  of  the  mentioned  characteristics: 
therefore,  a method  Co  measure  connectivity  was  devised.  Connectivity 
is  a measure  of  the  degreee  of  entanglement  through  petioles  and 
stolons.  Zero  connectivity  le  expected  In  macs  where  plants  ore  not 
touching  each  other,  and  the  plants  are  not  connected  by  stolons.  The 
principle  of  the  devised  connectivity  test  Is  chat  os  conneetivicy 
Increases,  more  wel^t/unlc  imt  area  Is 


required  to  submerge 


IS 


non-conneoced  plants  {Pa).  Sons  mats  had  similar  buoyancy  values  but 
dlflerent  connectivities.  To  ascertain  easily  the  affect  of 


Besolts  and  Discussion. 
Elam  Mass  Distribution  Histograms. 


standard  deviation 


the  nuaber  oF  samples  (Press  et  al.,  1966). 

Plant  size  vas  represented  by  six  classes  and  various  subclasses 
(Table  3.0):  this  classification  sysceo  smy  have  Co  be  expanded  by  Che 
addition  of  other  subclasses  Co  fit  ocher  population  groups  not 
sampled.  Class  0 represents  plants  chat  veigh  less  than  0.1  kg.  This 
class  was  subdivided  into  subclasses,  because  the  largest  and  scaliest 
plants  In  a given  population  fell  elthln  only  two  adjoining  classes 
(e.g.  Che  ssallesc  plant  fell  into  class  0 and  the  largest  In  1). 

Class  5 represents  large  plants  that  eeigh  over  0.7  l«-  The  class 
structure  follows  a linear  increase  in  nass  until  class  3.  Classes  3 
through  5 allov  for  greater  nass  differences  because  larger  and  older 
plants  are  structurally  different  from  youi^er  and  stialler  plants. 
Larger  plants  usually  have  more  leaves  and  fewer  stolons,  and  they 
have  part  of  their  nass  In  decaying  leaves  that  eventually  slough  off. 


Table  3.( 
Class 


Subc lass 


constructed  with  these  classes  were  normally 
distributed  (Figures  3.1.  - 3.3).  end  these  mss  distributions 


according  to  age  of  culture  or  habitat.  After  the  data  from  the  i 
unsyiBMtrlcsl  group  (Figure  3.4)  was  studied,  it  was  determined  tl 


be  represented  by  more  than  one  group  of  plants. 
IS  0 plants  and  the  others  by  higher  classes 
;roup  of  plants  was  not  well  established  because 
controlled  by  regular  harvests.  The  histogram 
' hyacinth  size  could  greatly  differ  from  one 
md  even  within  the  same  location  if  the 
lot  well  established. 


(Figure  3.S).  Th 

Illustrate  that  e 
location  to  onoth 
population  is  stl 
Rhizome  Length. 

Data  on  rhizome  length  from  all  locations 
graphed  together.  RhlztHBe  length  is  linearly  correlated  to  plant 
weight  (Figure  3.6).  Ungth  varied  from  1 to  21  centimeters.  Two 
least  square  equations  that  relate  rhizome  length  to  plant  lass  wen 
evident:  one  for  plants  that  showed  considerable  weevil  daioage  or  me 
concentrated  in  the  rhizomes  and  one,  for  healthy  plante  with  ai^le 


Biomass  density  In  this  study  varied  from  8 to  20  kg/m®  (wb).  ond 
population  density  varied  from  40  to  300  plants/n®.  However,  no 
apparent  relationship  exits  between  the  two  types  of  density  (Figure 
3.7).  One  estiimte  of  population  density  is  the  quotient  of  the 
biomass  density  (kg/m  ) of  a given  aggregation  by  its  average  plant 
mass  (average  of  7 plants).  This  simple  relation  functions  well  for 


suo|)eAJSsqo  io  iaqum)^ 


suoiiBAJSsqo  >0  jaquinrg 


suo|)eAJ3sqo  to  jsqiunN 


suofiBAiSsqo  10  JoqtunN 


'mBus|  suioziMd 


population  densities  loss  than  100  plants/n^,  but  It  is  very 
unreliable  for  oggregatlons  of  snail  plants  or  large  populations 
{Figure  3.8}.  Snail  plants  usually  Indicate  a population  that  Is  not 
veil  establishod,  in  which  case,  a reliable  value  for  plant  nass  can 
be  difficult  to  obtain  (Figure  3.4).  A more  accurate  prediction  of 
plant  number  vas  established  (Figure  3.9).  The  relation  was  developed 
from  Tucher's  (1981)  observations  that  bioiatss  density  Increased  with 
petiole  length.  The  least  square  relation  Is  a nonlinear  relation 
involving  the  logarithm  of  avorago  plant  height.  Biomass  density  and 
average  plant  length  are  the  independent  variables.  Both  of  these 

they  should  not  be  used  for  long-rooted  plant  systems. 

Connectivity  and  Buoyancy. 

Hat  buoyancy  in  this  study  varied  from  20  to  240  Pa  and 
connectivity  from  20  to  180  Pa  (Table  3.1).  Penfound  and  Earle 
(1948),  on  the  other  hand,  found  that  loading  capacity  varied  from  177 
to  430  Pa.  The  larger  value,  430,  prohcdily  represented  data  from  an 
extremely  compact  aggregation. 

Mat  buoyancy  did  not  correlate  well  with  standing  mass  or 
population  density,  although  buoyoncy  appeared  to  increase  somewhat 
with  areal  mass  density.  The  relation  between  mat  buoyancy  and  areal 

ratio  of  connectivity  to  mat  buoyancy  (Figure  3,10),  The  least  square 
relations  ore  presented  in  Table  3,2,  The  highest  buoyancy  at  any 
given  areal  mass  donaicy  occurred  when  connectivity  ratloa  were  less 
than  0.4,  and  the  lowest  rat  buoyancy  occurred  when  the  connectivity 


. Lu/s)ue|d  '/(jisusQ  uoijeindod  pajaipaid 


ed  'AsueXong 


connectivity  ratio  because  the  ratio  reflects  the  type  and  condition 
of  plants  cofoprlslng  a oat.  High  connectivity  ratio  Indicates  the 
plants  are  sui^rted  by  petiole  entanglement  or  are  Interconnected  by 
stolons,  while  a snail  ratio  indicates  many  petioles  are  floated  and 
Individual  plants  are  supporting  themselves.  Petioles  with  floats  and 
plants  of  low  connectivity  ere  emnnon  when  a oulturo  Is  periodically 
harvested,  because  open  conditions  lead  to  float  formation  In  plants 
of  any  size  (Rea,  1920).  Low  connectivity  results  whan  plants  are 
separated.  Floats  increase  nac  buoyancy,  because  the  specific  gravity 
of  a float  Is  0.136  or  ono-Flfth  the  density  of  a regular  petiole. 


Table  3.2  Relation  between  oat  buoyancy  (MB)  and  standing 
blosass  density  (BO) 

Equation  Connectivity  ratio  r 


MB  s -262.0  * m.5  Jn(  8D) 
MB  = -235.8  ♦ 139.2  ln(  BD) 
MB  V -212.3  ♦ 141.2  ln(  BD) 


Two  values  for  buoyancy  were  collected  for  lor^-rooted  plants, 
and  they  were  smaller  chan  those  expected  for  short-rooted  plants  with 
similar  connectivity  and  standing  densitiee.  This  is  expected  because 


long-rooted  plants  have  around  261  none  biomass  underneath  the  eater 
Suroerv- 

Vacer  hyacinths  consist  primarily  of  eater,  and,  although  comples 
structurally,  they  can  be  defined  by  loaves,  rhizomes,  roots  and 
stolons.  Floatless  petioles  help  support  the  plant  under  dense 
conditions  by  entanglement  with  leaves  from  other  stats.  Rhizome 
length  Is  correlated  to  plant  nass,  and  root  length  to  nutrient 
availability. 

distributed  plant  oass,  ehlle  aggregations  of  small  plants  or  plants 
subjected  to  regular  harvests  can  consist  of  many  types  of  plsnt 

Population  density  is  correlated  with  plant  length  and  standing 

Connectivity  and  mat  buoyancy  are  independent  properties  thnt 
define  the  mechanics  of  the  vertical  stability  of  nmits  and  plants. 


OUtTER  IV 

VISCOUS  AND  ELASTIC  PROPERTIES 


A aysten  is  a collection  of  inceractlre  eleoonts  for  «hich  there 
are  cause-and-effect  relationships  enong  the  varlahlos.  A eater 
hyacinth  system  can  be  subjected  to  friction  and  elastic  forces 

driving  forces  of  wind  or  towing  devices,  and  drag  from  the 
interaction  of  plants  with  water.  Therefore,  probable  system  elements 
are  related  to  body  forces  (inertial  and  gravl  talonnl)  and  surface 
forces  (pressure  and  shear).  To  represent  the  elastic  and  frictional 
forces  in  form  of  system  elements,  constitutive  equations  are 
required.  They  depend  only  on  the  characteristics  of  the  physical 
imterlBl  that  cofsprlses  an  element  and  on  its  geometry,  and  they 
describe  basic  relationships  between  two  variables  (e.g.  Hooke's  law 
in  solids  and  Newtons  viscosity  law  in  fluids). 

Hydrodynamic  Drag  Forces 

The  purpose  of  this  section  is  to  develop  the  constitutive 
relationships  between  drag  and  towing  velocity  os  a function  of  size 
of  an  aggregation  of  water  hyacinths  and  plant  size  and  shape. 


Review  of  I.i  lernciir*^. 


Kootfi,  rhlzonea  and  lower  petioles  of  water  hyacinths  are 
sutwnerged  la  water.  When  they  are  placed  In  flow,  they  are  subjected 
to  drag.  Drag  Is  composed  of  two  coeponencs:  friction  drag  resultlr^ 
from  tangential  stresses  caused  by  friction  of  the  fluid  against  the 
Interface,  and  pressure  or  forn  drag,  resulting  fron  pressure  forces 
exerted  by  the  fluid  nonial  to  the  surface  of  the  body.  The  latter  Is 
predominant  at  high  Reynolds  numbers.  Drag  on  a body  can  be 
determined  analytically  for  slaple  shapes  and  flows  and  numerically 
for  slightly  more  complicated  geemetries.  but  It  is  usually  determined 
experimentlally  for  bodies  in  turbulent  flow.  Custcxiairlly,  the 
formula  for  overall  drag  Is  written  in  the  following  form: 


where  D.  Is  drag  force.  A Is  a the  area  of  the  body  projected  In 
the  direction  of  flow,  e Is  density  of  the  fluid,  C,  is  a drag 
coefficient  and  V Is  the  velocity  of  the  fluid  relative  to  the  body 

A similar  drag  formula  was  Initially  proposed  by  Newton,  but  the 
original  proportionality  factor  (instead  of  C.}  did  not  always 
coincide  with  experimental  data.  This  problem  was  resolved  when  a 
function  of  Reynolds  number  was  used  to  replaced  the  original 
proportionality  factor  (Prandtl  and  Tletjens,  1934}.  Reynolds  number 
enabled  calculation  of  drag  in  a Newtonian  fluid  on  similar  bodies 

velocity.  The  relation  between  C.  and  Reynolds  number  must  be 


34 


y experluentially.  Furthermore, 


determined  for  turbulent  flow  oi 
experiments  proved  that  when  the  total  drag  on  a body  consists 
exclusively  of  pressure  drag  (l.e.  for  high  Reynolds  numbers),  the 
drag  coefficient  la  large  and  Independent  of  Reynolds  number.  In  the 
case  of  friction  drag,  the  drag  coefficient  Is  a function  of  Reynolds 
number,  and  It  Is  usually  small  In  value. 

Rhizomes  are  the  principal  projections  Into  water:  they  are 
cyllndrlcally  shaped,  thie  to  the  shape  of  rhizomes,  water  hyacinths 
In  flow  can  be  modelled  from  flow  analyses  for  a scries  of  cylinders 
(e.g.  In  a heal  exchanger).  Also,  Infonmttlon  concerning  the  flow 
properties  of  water  hyacinth  can  be  ascertained  from  a review  of 
studies  of  flow  retardance  by  aquatic  plants,  which  are  assumed  to 
represent  roughness  elements.  These  topics  will  be  reviewed  next. 

Flow  normal  to  a bonk  of  cylinders.  Since  tbe  shape  of  rhizomes 
can  be  assumed  to  be  cylindrical,  the  following  Infomation  on 
crossflow  over  a hank  of  cylinders  Is  pertinent  to  the  description  of 
water  hyacinth  behavior.  Crlmiscn  (1937)  correlated  shear  stress,  r 


encountered  by  flow  through  a bank,  and  F Is  a friction  factor  which 
Is  a function  of  Reynolds  number  and  the  ratio  of  transverse  to 
longitudinal  spaclngs. 

He  found  that  friction  coefficients  for  different  arrangements  of 
cylinders  were  functions  of  Reynolds  number  (based  on  cylinder 


restrictions  sncouAcered 


dloneter),  the  number  of  major 

The  number  of  major  restrictions  depended  on  Che  conflsuracion  of  Che 
cylinders:  for  in-line  arrangemencs,  p equaled  Che  number  of  rows; 

, with  oininum  flow  area  In  diagonal  incercyllndrlal 
spaces  or  openings  between  cylinders,  p equaled  Che  number  of  rows 
minus  one.  Generally,  the  friction  factor  was  constant  until 
Reynolds  nus^r  exceeded  10.000.  Above  10,000,  the  general 
characteristic  curve  for  friction  factor  Inclined  downeard. 

Later,  Che  equation  for  total  head  loss  due  to  friction  was 
refined  to  resemble  the  friction  factor  from  a Moody  diagram  (Equation 
4.4]  (l^ys  and  London,  19&4  and  HcQulston  and  Rerker,  1982): 

Total  head  loss  (Kays  and  London,  1964}  = 


r*v2 

» ■ ' A^  2g  • 

where  F is  Moody  friction  factor.  V is  mean  velocity,  _=  , 

A Is  total  friction  area,  A Is  total  flow  cross  sectional  area, 
and  if  Is  flow  length.  Both  A and  if  are  measured  from  the  leadlr« 
edge  of  Che  first  cylinder  row  to  the  trailing  edge  of  the  last  row  of 
cylinders,  and  R equals  the  hydraulic  radius.  Hydraulic  radius  for 
flow  normal  Co  cylinders  Is  cylinder  diameter. 


s often  presented  i 
>y  dlnenslonal  analysis  (Lanb,  1945): 


here  a is  a characteristic  length,  u Is  klnomtic  viscosity,  V U 
average  velocity  at  a cross  section  and  n Is  an  eicponent  obtained 
through  experlnentatlon.  Reynolds  number  is  defined  by  , For 

nsO,  the  fonrila  that  Is  often  used  in  hydraulics  for  the  case  of 
turbulent  flow  through  pipes  Is  obtained  It  Is  sinilar  to  Equation 
4,0.  Lanb  (1945)  suggested  that  the  value  of  ra  appears  to  depend  on 


Without  exception,  studios  of  flow  retordance  by  aquatic  plants 
have  assumed  that  plants  In  channel  flow  act  os  roughness  eletsents, 
Basic  ranges  of  Flow  reglises  for  roughness  elements  have  been 
established  on  the  basis  of  height  and  spacing  of  roughness  elenents 
(Schllchtlng,  1979),  These  are  1)  the  hydraulically  smooth  regime, 
where  the  size  of  the  roughness  Is  so  small  that  all  protrusions  are 
contained  within  the  viscous  sub-layer:  2)  the  transition  regime, 
where  protrusions  Influence  the  viscous  sub-layer,  and  the  resistance, 
compared  with  a smooth  pipe.  Is  due  mainly  Co  form  drag  experienced  by 
the  protrusions  In  the  turbulant  boundary  layer;  and  finally,  3)  the 
completely  rough  regime  where  all  protrusions  extend  outside  Che 
viscous  sub-layer  and  the  resistance  is  due  mainly  to  form  drag. 

These  regimes  are  determined  by  the  following  criteria: 


Hydraulically 


c.  Conplctely  rough. 

is  the  friction  velocity  defined  by  v*r  /p,  R Is  the  hydraulic 
radius;  j — Is  often  referred  to  as  the  oall  Reynolds  nunber,  and 


is  shear  stress  at  the  rail.  Friction  factor  {Oarcy-Welslach)  1 
related  to  shear  stress  by  = vf/B  V.  Furthermore,  It  tvis 

been  related  experimentally  to  Hanning's  n (Strlckor,  1929): 
n = CI/26)^  k,j6S) 


• 1th  units  n'^'^s  *,  and  la  roughness  In  meters  such  that  65X 

of  the  material  Is  of  lesser  size.  The  equation  is  valid  for  d < -^ 
<d000.  Hanning's  n Is  a roughness  coefficient;  Its  value  Increases 
with  Increasing  roughness. 

In  reference  to  flow  over  aquatic  plants,  Roe  and  Palmer  (1949) 
demonstrated  by  experiments  that  the  degree  to  which  bottom  channel 
vegetation  retards  water  flow  depends  largely  upon  the  degree  to  which 
vegetation  Is  bent  and  flattened  by  the  flow,  which  In  turn,  depends 
mainly  upon  physical  characteristics  of  the  vegetation,  Its  iianner  of 
growth,  and  the  velocity  and  depth  of  flow.  The  plants  were  different 


sufficiently  deep  and  fost  to  bend  over  and  submerge  the  channel 


vegftcaclon,  & greater  velocity  resulted  duo  to  a greater  slope  ol  the 
energy  grade  line  resulted  in  loser  resistance  to  rioo.  Hanning's  n 
values  for  short  and  ncdlun-helght  vegetation  appeared  to  decrease  to 

to  cover  the  planes.  Values  of  a for  different  tall  plane  species 
varied  widely  for  all  depth  of  flows.  In  regard  to  population  density 
for  a specified  flow  depth,  a dense  cover  of  vegetation  (300 
plents/a^)  offered  twice  as  much  resistance  to  flow  as  a less  dense 
cover  (17  plants/m  ),  The  authors  developed  a graphical  method  for 
the  determination  of  channel  cross  section.  To  do  so,  they  plotted 
Hannings  n values  versus  VR  (velocity  times  hydraulic  radius)  because 
they  perceived  resistance  to  be  a function  of  the  degree  of  flattening 
of  the  vegetation,  which  is  influenced  by  velocity  and  depth  of  flow. 
As  a result,  they  found  chat  until  a given  vegetative  channel  llnli% 
became  well  flattened.  Che  value  of  VR  indicated  Che  degree  of 
resistance  the  lining  offered.  However,  when  the  vegetation  was 
prone,  or  nearly  so,  and  well  subsicrged,  n became  practically  constant 


On  the  ocher  hand,  Kouven  and  Unny  (l€ff3)  dotemlned  two  values 

regimes.  They  found  chat  friction  for  the  erect  and  waving  regime  was 

roughness,  friction  factor  was  a Function  of  VR.  The  experimental 
reeulcs  were  in  the  form  of  graphs.  This  apparent  discrepancy  in  the 
two  investigations  could  relate  to  the  amount  the  vegetation  was 
submerged  or  to  the  depth  of  flow.  Although,  It  appears  that  Kouwen 
and  Urmy  were  correct  In  their  conclusion  because  1)  flow  theory  of 


a fully  rough  regime  In  a channel  predicts  the  drag  coefficient  to  be 
Independent  of  Reynolds  number,  so  drag  varies  with  and  is 
independent  of  viscosity:  and  2)  erect  vegetation  is  more  likely  to 
cause  a fully  rough  flow  regime  because  the  roughness  elements  are 
bigger  to  provoke  larger  values  for  the  wall  Reynolds  number. 

ICany  Investigators  have  found  that  Mannings  n or  friction 
coefficient,  increases  with  the  amount  of  aquatic  biofass  present 
(Bogart,  19d9:  Kouwen  et  a].,  1969;  Powell,  1978:  Petryk  and  Bosoajran 
111.  1975:  Shlh  and  Rahi,  1982  and  Stephens  el  al..  1563).  Petryk  and 
Bosmjran  III  (1976)  defined  density  of  aquatic  vegetation  by  the 
vegetative  area  per  unit  length  of  channel  per  unit  area  of  Mow  as 
aitpressed  by  the  following  equation: 

t A, 

Density  of  vegetation  = 


where  Aj  = projected  cross  sectional 
direction  in  a channel  of  ler^th  aj 
The  effect  of  plant  density  on  n< 
roughness  spocings.  Powell  (1946)  stui 

roughness  height  was  an  Important 
roughness,  and  it  was  constant  regardless  of  the  value  of  the 
roughness  width  to  height  ratio.  Sayre  and  Albertson  (1961)  sti 

il  analysis  to  develop  the  following  spacing  parameter: 


long!  tudliusl  spacing  of 


Spacing  parameter  - ^ : (d-5) 

! X - longitudinal  spacing  of  roughness  teffle:  h 


roughness  width  and 


spacing  lei^th. 


The  flow  was  considered  Co 


IV  near  a rough  boundary. 


the  ooi^lecely  rough  range,  Ebca 
rman-Prandtl  concepce  of  turbulent 
the  following  equation  was  fitted  to 


whore  C Is  the  Oiezy  discharge  coefficient  = R*'^/n  , y Is 

normal  depth,  Is  an  ewperinental  constant  dependent  on  roughness 

spacing,  h Is  height  of  roughness  elemenc,  and  g is  acceleration 
of  gravity.  Similarly,  using  the  equation  for  logarithmic  velocity 
profile  near  a turbulent  wall,  Kouwen  et  al , (ld69)  obtained  this 
equation  for  flow  velocity  In  a vegetated  channel: 


where  Cj  and  Cj  vary  with  the  vegetative  density  and  the 
flexibility  of  the  vegetative,  A is  the  gross  cross  sectional  area 


of  Che  channel  and  A la  Che  area  of  vegetative  cross  section 

for  application  to  water  hyacinths  are: 

a.  Theory  and  experiments  have  shown  that  bigger  plants  aj 
denser  vegetation  produce  a larger  drag  coefficient. 


Che  bonk,  specifically,  cylinder  length,  bank  length  and  bank  width 
Influence  friction  facCora. 

c.  For  Reynolds  nunber  (based  on  tube  dlanecer)  less  chan  10, 0i 
the  friction  factor  of  a bank  of  cylinders  Id  oross  flow  is 
Independent  of  Reynolds  nus^r. 

d.  Transverse  and  longlcudial  spaclnga  of  cylinders  In  cross 
flow  and  of  roughness  projections  in  flow  affect  drag  coefficient. 

e.  Plants  chat  are  flatten  In  flow  offer  less  reslscence  to  f 


n vegetation  ii 


drag  coefficient  Is  Independent  of  Reynolds  number,  and  it  Is 
dependent  on  a relative  roughness  parometer. 

In  this  section,  an  ewpressicn  for  drag  valid  for  a rectai^lar 
aggregation  of  plants  of  any  population  density  or  spacing  and 
frictional  area  In  flow  Is  developed  In  this  section.  Neither 
nunerlcal  nor  analytical  solution  appeared  feasible  for  the  problen  of 
water  hyacinths  in  steady  flow  due  to  the  assuoptlons  of  1)  rhizomes 
and  roots  cause  turbulence,  and  2)  the  mat  deforms  when  moved. 
Therefore,  an  analytical  solutloo  based  on  rigid,  rough  flat  plote  Is 
not  appropriate.  Consequently,  a drag  ns 
Newton's  drag  equation  (Equation  4.0},  ai 

To  begin,  a model  for  drag  must  be  dimensionally  correct.  The 
drag  equation.  Equation  4,0,  was  used  after  it  was  rodlfled  by  Che 
Incorporation  of  Che  number  2 into  Che  drag  coefficient. 


n rigid,  rough  flat 
del  was  developed  be 
id  drag  coefficients 


Df  = p * 

« equals  the  drag  force  (Newtons), 

:r  approxlDOtely  equal  to  1000  . A 


(4.6) 

equals  the  projected 


area  In  flow  (d  ),  C.  equals  a nondlisensfonal  drag  coefficient,  an< 
V Is  the  towing  velocity  ( a/s).  Although  Equation  4.0  Inplies  C 
is  a function  of  Reynolds  nunber,  evidence  In  the  literature 
Indicates  that  the  drag  coefficient  for  an  erect,  rigid  projection 
like  a rhizoae  of  water  hyacinth,  is  Independent  of  Reynolds  nunber. 
If  Indeed,  C,  of  water  hyacinths  is  Independent  of  Reynolds  number 
the  towing  velocity  would  enter  the  owidel  only  as  V^.  To 
substantiate  this,  the  Reynolds  number  was  calculated.  Dlamotors  of 
rhizomes  range  from  0.01  to  0.03  m.  The  Reynolds  number  for  water 

Is  calculated  below: 


However,  for  0.03  m diameter,  the  Reynolds  nunber  Is  approximately 

have  Reynolds  numbers  below  the  lower  limit  at  which  drag  of  cylinders 
in  cross  flow  is  dependent  on  Reynolds  number. 

To  use  the  drag  equation,  the  friction  area  must  be  defined.  The 
total  friction  area  was  assumed  to  approximately  equal  the  projected 


nuDib«r  of  planes  (N  } 


«>  of  a rhizooe  (A  ) oulcipHed  b; 
an  aggregation  or 


D equals  Eha  average  diameter  of  Che  rhizomes  and  L , average 

coefficient  since  drag  coefficients  are  only  similar  for  geometrically 
similar  bodies  of  the  sane  orientation  In  relation  to  the  free  stream 
velocity.  Therefore,  a nat  which  Is  wider  or  has  more  plants  along 
the  flow  direction  will  have  a different  drag  coefficient  than  a 
narrower  mat.  The  shape  factor  enters  the  cylindorieal  crossflow 
models  as  p (Equation  4.4)  and  nondlmenslonally  as  A/A  (Equation 
4.5).  For  a rectangular  mat  of  water  hyacinth,  the  simplest 

width.  W.  to  length,  L. 

For  this  study  the  transverse  and  longltudial  spaclngs  of  water 
hyacinths  within  a mt  was  assumed  uniform.  As  a result,  the  drag 
coefficient  was  considered  to  be  a function  of  a shape  factor  and, 
possibly,  Reynolds  number.  Then, 
hyacinth  becomes. 


(4.8) 


'HilB  proposed  equation  li 
of  eater  hyacinch  without  the 

Materials  and  Methods. 


o 1) 


Che  hypothesis  that  the  drag  coeFficlent  of  a water  hyacinth  sat  is 
Independent  of  Reynolds  number  and  2)  decerolne  tho  empirical 
relation  of  to  mat  geometry. 

Aggregations  of  water  hyacinths  composed  priisarily  of  small, 
medium  or  large-sized  plants,  were  cowed.  Small  plants  averaged  25  cm 
In  length,  medium  plants  averaged  61  cm,  and  large  plants  averaged  90 

on  the  mat  before  it  was  moved  to  the  testing  area,  ffater  hyacinths 
were  obtained  from  the  University  of  Florida  Swlno  Research  Unit. 
Blven's  Arm.  and  Lake  Alice,  The  lakes  are  located  near  the 
University  of  Florida.  After  cospletlcn  of  each  experiment,  seven 
plants  were  measured  for  irmss,  rhizome  length  and  diameter,  largeat 
petiole  length,  plant  length  and  root  ler^ch.  Root  length  varied 
according  to  habitat  and  reused  from  2 to  24  cm.  The  presence  of 
stolons  and  floats  waa  noted. 

Aggregations  of  plants  from  the  Swine  Unit  were  relocated  to  the 
tasting  area  in  this  manner.  First,  a 2.44  by  1.22  m area  of  plants 
was  isolated:  then,  a frame  of  wood  and  fence  wire  was  positioned 
underneath  the  plants.  Lifting  the  frame  removed  the  plants  from  the 


positions.  Because  Che 


rrantr-llke  stniccure;  a po»ercd  trolley;  rater  velocity 
nd  instruMBtatlon  (Figure  4.0). 


bridges. 

rentally 


‘event  plants  Fron  cllinbing  ovei 


volume  (Figure  4.1). 


Before  some  of  the  tests,  a grid  of  thin  nylon  t»tne  ras  placed 
over  the  eat  of  plants  for  subsequent  observation  of  plant  rovement  in 


suspended  fron  the 


cantilever  heans  which  ware 
e and  attached  to  the  lower 
frame  contained  plants.  Towins  velocity 


To*lng  speed  was  increased  froie  0-1  to  0.6  m/s  until  the  plants 
started  to  roll  under  the  conialment  structure.  When  this  happened, 
Che  test  was  ended,  and  the  sample  was  subdivided.  Some  of  Che  plants 
were  removed  to  obtain  another  mot  size.  Retooved  plants  were  drained 
and  weighed  to  determine  the  quantity  of  biomass  in  the  sample. 

Towing  tests  were  repeated  on  sequentially  ssmller  netCs  until  all  mt 


Force  and  velocity  data  for  each  test  were  regressed  acoordir^  Co 

Df  = (4.9) 

where  equals  total  drag  force  in  Newtons.  C is  a constant  of 

proportionality  equal  to  the  product  of  [p  N D L C.)  from  Equation 
4.8  and  V equals  cowing  velocity.  Drag  force  was  determined  to  be 
the  force  on  Che  cantilever  beasEs  after  the  plants  had  stopped 
accelerating  and  the  forces  due  to  the  acceleration  stage  of  the  test 
had  diminished  to  a slight  fluctuation.  The  parts  of  Che  cowing  tests 
were  labelled  the  unsteady  and  steady  stages  (Figures  4.2  and  4.3). 

The  data  were  analyzed  for  error.  Random  or  experimental  error, 
can  be  quantified  when  its  distribution  ai^rooches  a normal 
distribution  (Herschy.  1978).  The  principal  source  of  experimental 
error  was  the  data  logger.  Drift  of  the  sec  point  of  the  data  logger 
was  i 0.083  mv,  which  translates  into  an  error  of  h 0.9  N.  Skewness 
of  the  force  data  during  the  steady  state  conditions  of  a towing  test 


N‘»3J0d 


N 'aojoj 


5 calculac«d  for  t* 


the  distribution  Is  nomal  i accordingly,  Che  force  distribution  can  be 
considered  Co  be  noml.  A hlscosram  of  the  drag  force  distribution 
1s  shown  In  Figure  4.4.  The  dispersion  of  drag  force  about  the  tlieo 
average  was  tncasured  In  standard  deviation  since  the  distribution  was 
noroBl , Standard  deviation  of  drag  force,  D^,  was  calculated  using 
data  froit  the  two  tests;  while,  standard  errors  of  (Equation 

coefficient,  C,  (Equation  4.11)  were  computed  for  every  test.  Since 
the  Intematlonal  Standards  Organisation  has  roconnended  that  the 
probability  level  used  should  be  at  the  95  percent  level  (Herschy. 
1978),  error  was  reported  using  tho  corresponding  t statistic 
multiplying  either  Che  standard  error  or  deviation. 


Standard  Error  of  Estimte,  S (Herschy,  1976), 


(4-10) 


where  d is  the  deviation  of  the  actual  value  from  the  computed  value 


Standard  Error  of  the  Cbefflclenc  (Herschy,  1978)  Is, 


- (I  v2)2- 


(4.11) 


siuaiuajnses^  iaqiunN 


Spurious  errors  were  discarded  {Flsure  4.3)  when  data  lay  more 
depends  upon  the  sample  population  size  (Benedict,  19S5).  For 

Systematic  error  la  the  flow  tests  could  be  related  to  the  deslgsi 
of  the  teat  frame.  Tines  retained  plants  within  the  test  frame.  The 
resistance  to  flow  that  was  considered  small,  because  the  diameter  and 
number  of  tines  were  small  relative  to  those  of  the  rhizomes. 


velocity  of  0.1  m/s  and  3.0  N at  0.5  n/s.  The  drag  on  the  saaller 

required  In  any  towing  system,  so  their  exclusion  from  the  experiments 
was  not  deemed  necessary.  Also,  force  measuramencs  could  be 
eusceptlble  to  systemic  error  If  the  eallbraticm  of  the  strain  gages 
was  loaccurate.  The  gages  were  calibrated  twice  using  different  input 
voltages,  spring  scales  and  ranges  of  force,  and  the  calibration 
factors  reimined  the  same.  Finally,  another  source  of  systematic 
error  ie  the  number  of  significant  digits  the  data  loggsr  uses  to 
score  force  readings.  Too  few  and  very  small  forces  would  not  be 
recorded  accurately  by  the  data  logger. 

Results  and  Discussion. 

Experimental  observations. 

Ac  the  beginning  of  a towii^  test,  plants  tended  to  bunch 
together,  causing  the  length  of  Che  sac  to  decrease.  Analysis  of  the 
movement  of  grid  lines  on  top  of  the  hyacinths  demonstrated  mechanisms 
of  compaction.  When  a test  began,  the  grid  lines  at  Che  leading  edge 


of  Che  nac  reealned  stationary  chose  at  the  rear  moved.  Once  Che 
whole  mac  reached  constant  speed,  Che  grid  lines  were  equally  spaced 

Compaction  Increased  with  toeing  velocity  until  "roll  under". 

point.  Che  plants  on  the  leading  edge  started  to  roll  forward  under 
Che  water,  ftice  under  water,  the  plants  tended  Co  flatten  out 
underneath  Che  other  plance  aoving  on  Che  surface.  In  this  way,  Che 

rolled  under  when  they  were  positioned  at  the  leading  edge  (Figures 
d.5  -4.7).  Ac  "roll-under",  the  force  on  the  cantilevers  greatly 
Increased,  and  the  portion  of  the  force  - velocity  graph  pertaining  to 
"steady"  was  absent  (Figures  4.8  and  4.9). 

results.  Data  from  repeated  tests  fell  on  Che  sudel  line  (Figure 
4.9):  hence.  Che  testa  could  be  duplicated  under  similar  conditions. 

Vhen  plants  are  stationary,  their  roots  hang  downward  In  the 
water;  yet.  they  are  very  flexible.  Underwater  observations  showed 
the  roots  to  be  completely  prone  at  a towing  velocity  of  0.17  m/s  and 
waving  hack  and  forth  at  lower  speeds.  Once  the  roots  wsre  prone,  the 
component  of  drag  was  due  to  rhizomes. 

Results  of  reeresslon  analysis. 

A model  fits  experimental  data  provided  the  standard  deviation  of 
Che  dependent  variable  fells  in  the  range  of  the  computed  standard 
error  of  the  predicted  value.  To  determine  the  efficacy  of  the 
proposed  model,  time  average  values  of  the  dreg  force  were  plotted 
egalnsc  the  square  of  the  towing  velocity  for  two  sample  testa.  The 


When  nace  were  loved,  they  compacted,  amount  of  compaction 
increaaed  with  cowins  volocity.  At  0.1  m/a,  compaction 


Fljure  1|.6  The  nec  wae  coved  ec  0.3  e/i.  end  it  conpeeced  appcoxuacel; 


Figure  4.7  The  net  MS  novlng  sC  0.4  a/a,  end  1C  co»i«eced  approxlnacely 
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(shown  as  S£) 


plotted 


deviations  fell  in  the  range  of  the  standard  error  (.9S  confidence 
level).  'The  fact  that  the  ranges  of  the  standard  deviations  of 
at  "roll  under  " velocity  did  not  fall  between  the  ranges  of  its 
standard  error  shows  that  the  flow  behavior  of  the  water  hyacinth 


2 .(Equation  4.9)  were 


very  quickly;  In  those  cases,  few  date  points  were  obtained.  Morever. 
82  percent  of  the  C values  had  less  than  ten  percent  relative  error. 
Error  was  not  consistent  In  each  eyperlnent  nor  between  experlnents. 

physclal  variations  In  the  nat  and  Irtdlvldual  plant  s 


Provided  that  Equations  4.6  and  4.9  accurately  describe  drag  for 
any  choice  of  plant  and  net  size  and  that  C.  is  not  a function  of 


N ‘aoioj  6ejQ 
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The  values  of  obtolned  by  this  expression  from  C are 

tabulated  in  Table  4,1  aloi%  with  values  of  the  shape  factor  and  othor 
variables.  Varlobles  and  were  calculated  from  expressions 
developed  eisplrlclaly  In  Chapter  3.  RhlzOTe  diameter  Is  difficult  to 

to  vary  little  within  a population  or  by  os  little  as  0.6  cn.  the 
following  classification  schone  was  used:  short  plants  have  1 cn 
diameters,  medium  plants,  2.  and  call  plants,  3 cm. 

To  assess  the  accuracy  of  the  experimental  drag  coefficient, 
experimental  error  of  Equation  4.15  was  calculated  by  the  following 
equation  (4.16)  (Horsehy,  1978); 


denotes  the  percent  ei 


le  drag  coefficient 

rage  Is  apprcxlsately  10 
h the  quality  of  water 


Density  of  the  water,  p . could  < 
(dirty  or  clean)  and  temperature;  1 


were  assumed  to  be  10*  and  3X  respectlvlty  since  chat  Is  approxlsntely 
the  error  Introduced  by  the  esplrlcal  expressions  of  Chapter  3.  The 
error  due  to  was  estimated  to  be  3*.  After  substitution  of 
those  values  of  percent  error  into  the  error  equation,  the  percent 
experlraenlnl  error  for  the  drag  coefficients  was  15*. 


deteroloed  to  be  nonlinobr  (Figure  4.12},  and  It  is  written  below. 


cantilever  beajss.  The  s> 
voltage  to  the  strain  ga« 


SD  the  etain  gagee  positioned  on  the 

could  not  be  explained  until  alter  videotapes  ol  the  experiments  were 
viewed.  The  plants  were  seen  bending  In  the  flow.  The  bending  motion 
favored  friction  drag  over  pressure  drag;  therefore  the  drag 


expected  (Schllchtlng,  1979) 


73 
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The  regression  equation  relating  C.  to  V/L  was  analyzed  (or 
experimental  and  dispersion  errors.  Dispersion  of  the  data  nround  the 
regression  curve  Is  given  by  tS  where  t la  the  t statistic  (.95 
confidence  level]  and  S Is  the  standard  error  due  to  the 
regression. If  the  regression  equation  Is  valid,  the  confidence 
Intervals  of  for  every  shape  factor  should  fall  within  the 

range  of  C.  given  by  the  standard  error  of  the  regression.  This  was 
true  except  for  W/L  equal  2 {Figure  d-I3].  However,  the  discrepancy 


-Hydrodynamic  drag  coefficient  of  water  hyacinths  is  Independent  of 
Reynolds  nuoiber  but  dependent  on  the  shape  of  a mt. 

-The  water  hyacinth  drag  equation  was  based  on  Newton's  drag  formula, 

rhizomes  and  the  square  of  velocity. 

Elastic  Characteristics. 

Elastic  defornatlon  of  water  hyacinths  occurs  when  an  external 
force  applied  on  a mat  compresses  plants  against  each  other.  I^rts  of 

rhizomes.  The  objective  of  this  section  Is  to  determine 
experimentally  the  elastic  modulus  nid  Poisson's  ratio  of  the  rhizome. 
Literature  Review. 

material  of  two  or  more  components  having  various,  usually  contrastlre 
physical  and  mechanical  properties.  Rater  hyacinth  rhizomes,  like 


)U9>9!|903  6ejQ 


>.  1961). 


have  ordered  acruccure  (Tarnopolskll  and  Kinds 
Hacro-structurally.  the  rhlzoiae  coneiats  of  four  regions:  i 


Clastic  projiertlea  are  aasuiaed  to  vary  principally  in  the  radial 

structure  of  rhizoaies.  The  elastic  properties  of  a corn  cob  are 
sinilarly  characterized  oi>d  have  been  extensively  studied  (Anazodo, 
1663;  Anazodo  and  Chikwendu,  1983:  and  Dilkwendu  and  Anazodo.  1961). 
The  methods  and  equations  those  investigators  used  were  adopted  in 
this  work  to  detemine  the  elastic  response  of  water  hyacinths  to 

equations  for  the  deteniinatlon  of  elastic  modulus  E and 
Poisson's  ratio  -p  of  rhlzones  under  quasi-static  radial  compression 
between  two  steel  flat  plates  using  a large  deformation  approxissttion. 
are  (Oilkwendu  and  Anazodo.  1961): 

£ _ 2 Pfl-  0^1  t2/Y^  - 6/g  <■  67  Y^/288  -175  Y^/211921  jgj 
r (1  -Y*/  8 e Y'’/192) 

where  P is  the  known  experliecntally  applied  force  per  unit  length  of 


of  the  cylinder; 


w Is  Che  nonml  epproach  of  the  flat  plate  to  the  center  of  the 
cylinder,  or  half  of  the  total  nortaal  dcformtlon  of  the  cylinder: 

X = In(VY)  -0.6; 

Y Is  Che  half-angle  of  contact  between  the  cylinder  and  Che  flat  plate 
which  can  be  calculated  by  the  followli^  eouatlon: 

, flnlVYl-  0.5Ul-Y^/a  ♦y'*/1921  *(5  1^/48-267  TVliaaOl 
(2/Y^  - 5/6  ♦ 67  Y*/288  - 476  7^24192) 

where  r - radius  of  the  cylinder; 


^ ^ Y^  ( 1-  Y^/6  > 

* 4(1  - y®/8  ♦ Y^/192) 
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Values  from  these  equations  which  ore  based  on  a large 

elastic  contact  theory,  and  It  was  concluded  that  Che  large 
deforaatlon  theory  was  the  aost  accurate  but  only  by  one  percent. 

When  two  cylindrical  bodies,  with  their  axes  parallel,  are 
pressed  in  contact  by  force,  P , per  unit  length,  they  nsike  contact 
over  a long  atrip  of  width  2d  lying  parallel  to  the  axis.  To 
determine  the  value  of  W . and  total  diametrical  compression.  Hertz 

are  as  follows  (Johnson,  198S) ; I]  The  surfaces  a 


2)  The  strains  are  small:  (f  <<  r.  3}  Each  solid  can  be  considered 


as  SJ1  elastic  half-specs:  d « r,  rf<<  length  of  the  cylinder;  i 
the  surfaces  are  f rlctlcmless. 

Johnson  {1565}  states  that  the  Hertz  theory  is  valid  for  sons 


slnllarly  to  that  li 


theory  to  predict  ci 


‘ d (CMksendu  and  Anozodo,  1984);  therefore. 

Finally.  aithoueJi  the  rhlzone  is  covered  with  roots, 
IS  friction  between  contacting  rhizoncs.  Hertz  theory 

only  if  the  elastic  constants  of  the  two  sBterlals  are  different 
{Johnson,  1985).  Therefore,  the  Hertz  contact  theory  can  be  used  for 

of  a rhizoDie  dionecer  coapreased  between  two  other  hyacinths  through 
the  Bid-points  of  the  contact  areas  was  used: 


flontacting  cylinders  or 


E*  = 


Hpterinls  and  Wathoda. 

Jn  this  section  the  procedures  used  Co  decernlne  Young's  Ten 
and  Poisson's  ratio  for  eoter  hyacinth  rhizones  ore  discussed, 
values  will  be  used  whenever  the  total  diametrical  corspresslon  o 
rhizoBie  Is  needed  (Equation  4.20). 


Rhizomes  wi 
Long  rhlzoBies  w 


■e  collected  from  Blven's  Arm.  a lake  where  water 
harvested;  consequently,  the  rhizomes  were  long, 
e collected  so  Chat  long  test  samples  could  be  cut. 

dimensions  are  typically  prepared  for  a compression  test  (Sitkei. 
1936).  Roots  were  left  on  the  rhizomes  CO  correspond  Co  an  actual 
compression  condition.  After  the  samples  were  cut,  rhizomes  were 

probably  represented  internodol  material.  The  Intemodal  area 
produces  the  roots  and  leaves.  The  amount  of  spongy  nmterlal  varied 
from  sample  to  sasqjle.  The  samples  were  subjected  Co  radial 
ooBTpression  on  an  Inscron  universal  cestlte  machlno  until  the  rhizome 
cracked.  Loading  rates  of  25.  SO  and  127  cm/min  were  used  to  check 
viscoelastic  tendencies. 

Typical  force-deforoatlon  behaviour  of  a radially  compressed 
rhlzcsse  eoCBposite  eshlbits  two  distinct  stages.  During  the  beginnlr^ 


Che  rblzooe  beneath  Che  spongy  efiCerlBl.  Unlike  nany  biological 
iBCerlala,  rhlzoDos  do  not  have  a blo-yleld  point:  although,  a 
definite  rupture  point  was  evident. 

Young's  Toodulus  was  calculated  fron  Equation  d.l9  with  values  of 

curve.  Poisson's  ratio  was  caiculated  by  Che  raethod  of  Chlkwendu  and 
Anazodo  (19S4)  and  the  experimental  setup  Is  briefly  reviewed  below. 
Successive  photographs  of  the  cross-section  of  each  rhizome 


by  Che  use  of  Equation  4.20. 

HesuUs  and  Discussion. 

Poisson's  ratio  was  successfully  determined  for  two  samples 

0.49  and  0.45.  They  are  considered  typical  for  a saterial  Chat 
approaches  chat  of  a liquid  {Nohsenln,  1978).  The  rhizomes  are  90 
percent  water;  therefore,  chose  values  are  expected.  Hence,  Che 
average  value,  0.47.  was  used  In  the  calculations  of  Young's  modulus. 

The  Young's  modulus  ranged  frost  7.6  to  11  with  an  average  of 
9.9  HPe  (Table  4.2).  Loading  rate  appeared  to  have  no  effect  cm  these 
values.  The  magnitude  of  the  moduius  was  between  chose  for  corn  cob 


{18  NPa) 


(1.4  HPa)  (Oilkwendu 


Loedlng  v«locicy 


a HYACINTHS  TOWS 


This  chapter  Is  concerned  vlth  developmenc  of  a oiachefiiaclcal 

that  eleiDcnts  of  the  syscen  and  their  interconnections  obey. 
Principal  ayscee  variables  are  eat  horizontal  velocity  and 
dlsplacofsent.  lIsinB  toeing  velocity  as  the  Input,  the  niodel  vas 
solved  nunerlcally,  and  calculated  values  of  towlr^  force  and  sat 
cosipaction  were  compared  to  expcrlsantal  data.  Further,  In  this 
chapter  the  physical  condition  of  the  system  prior  to  vertical 
displacement  or  "roll  under"  Is  described. 

Conceptual  Hodel, 

To  construct  an  accurate  nathesallcal  s 
their  Interconnections,  knowledge  of 
required.  At  the  beginning  of  this  study,  physical  m 
water  hyacinth  mat  and  their  interconnections  were  obscure  and  ill 
defined.  So  the  pattern  of  different  coisponent  interconnections 
sanlpulated  until  the  resulting  nathesatical  model  generated  aocum 
Input-output  relationships  for  different  boundary  and  input 

that  consists  of  a water  hyacinth  met  of  length  L,  and  width  V, 
a device  that  encloses  them  Is  shown  in  Figure  5.0.  The  distance 


Mechanical  System 


length  L of  the 


The  enclosing  device  is  analogous  to  the  test  equlpnent  used  in 
the  towing  eicperlnants  described  in  Chapter  4.0.  For  other  towli« 
devices,  different  kinds  and  arrangenent  of  eioments  could  cnodcl  the 
enclosure.  It  Is  represented  by  a spring  and  a mss  elODienc.  The 
spring  (representing  the  cantilever  beasts  on  the  toning  test 


elesienis  and  a 


i by  teo  types  of  friction 
i friction  eleatents  la 
water  hyacinths  moving  through 

ie  viscous  friction  element  is 


associated  with  sliding  friction  of  contacting  petioles. 

Initially  the  system  is  at  rest.  At  tine  t > 0.  th 
moves  at  velocity  Vg{t).  As  the  eat  begins  to  move,  it 
Vhen  compaction  ends,  the  entire  mt  moves  uniformly  at  t 
the  enclosure  device. 


Hass  elements.  To  represent  part  of  a physical  system  by  an 
ideal  translational  mass,  it  Is  necessary  that  all  of  the  r«rtlcles  be 
connected  together  and  chat  they  move  with  substantially  eoual  or 
proportional  valocicles  and  accelerations  (Shearer  et  al..  1971).  The 
constitutive  equation  for  an  ideal  mass  element  is  devlved  from 
Heeton's  second  law,  and  Inertial  force  equals 


(5.0) 


1 ddnoccs  node  number. 

IBS  of  the  enclosure  - H . 

.ss  of  water  hyacinth  raot  plus  nt 


Vj  = Velocity  a water  hyacinth  mat  at  node  1 with  respect  to 

gives  the  additional  force  required  to  accelerate  the  sass  of  fluid 

water.  The  additional  amount  of  water  issss  depends  on  the  shape  of 
the  moving  object  (Gerhart  and  Cross.  1985).  Water  hyacinth  rhizomes, 
the  principal  objects  in  flow,  are  similar  to  circular  cylinders.  For 

mass  due  to  water  equals  the  mass  of  fluid  displaced  by  the  circular 
cylinders  (Prandtl  and  Tletjens.  1934). 

element,  provided  that  an  algebraic  relationship  exists  between  the 
the  towing  experlmente  can  be  Idealized  as  one  stiffness  element.  For 


= K Xo. 


deflection  that  Is  given  by: 


= Xo  - X,, 


dlsplaceiBent  of  a water  hyacinth  aal  at  node  i.  Also, 
deflection  Is  defined  to  be  zero,  with  Xe  - X,  s 0. 
applied  to  the  cantilever. 

For  SBiall  deflections,  6.  of  a cantilever  beaio,  K 
calculated  as  follows: 


CS.2) 


Vhere  P = concentrated  load  at  free  end  of  beam, 
E - modulus  of  elasticity. 


h B beam  thickness. 

For  b B 5 cm,  h = 0.3175  cm,  t = 0.30  m , and  E {for  steel)  = 
four  cantilevers  beams  actlr^  conjunctively  so  the  theoretical  value 


Another  useful 


(8.3) 


;e  depends  oi 


The  period  o 

is  Independent  of  sagnl  tude  of  osclUaclons. 
equal  to  12.0  kg  and  K equal  to  12.000  n/a,  the  systea  natural 
frequency  Is  5 hz. 

Viscous  eleaent.  Drag  friction  Is  primarily  due  to  the  presence 
of  water  hyacinth  rhizomes.  If.  In  addition,  the  enclosure  device  1; 


equipped  with  objects  that  protrude  Into  the  wat. 
visclous  eleaents  may  have  to  be  added  to  Che  mo 
the  description  of  drag  due  to  rhlzotses  was  foun 
function  of  mat  velocity.  In  the  conceptual  modi 
the  velocity  of  node  1.  Therefore,  the  steady  f 
viscous  element  is 


, different  types  of 
to  be  a nonlinear 


Dp  = Crf  ''l 

The  absolute  value  symbol  is  required  to  change  direction  of  force 
with  change  In  direction  of  velocity.  Drag  coefficient  for  steady 
flow  was  found  to  bo;  = 0.-1-136  eiip(  0.906  -f-  ) (Equation  4-18). 

Drag  on  a mat  in  unsteady  flow  was  not  determined  experimentally, 
because,  several  physical  components  in  the  system  were  simultaneously 
active  (a.g.  Inertia,  drag,  and  contact  friction  elements),  and  drag 
coefficient  la  not  known.  Vfhen  a bluff  body  starting  from  rest  moves 
through  a fluid,  wakes  and  vertices  form  and  change  with  time. 

resulting  from  the  lengthening  of  these  wakes  (Panton,  1964).  Also, 
the  drag  coefficient  of  a group  of 


« runoclons  oJ 


with  acceleration  in  a turbulent  fluid  «as  found  ti 
particle  Reynolds  nunber  and  degree  of  turbulence,  but  not  of  the 
acceleratloo  (Klker  ejid  Ross,  1966  and  Torobln  and  Gauvln,  1961), 

Equation  5,4  while  a nu  Is  occeleratli^,  because  the  towlr^ 
experlnents  and  theory  described  In  the  previous  chapter  give  a 

Increase  of  wakes  Is  constant,  or  when  the  effects  of  systen 
acceleration  are  absent.  Therefore,  an  expression  for  C,  that  Is 
valid  for  unsteady  flow  conditions  was  obtained  through  nodel 
simulation. 

Contact  friction.  Another  friction  force  arises  because  plants 

condition  chat  occurrs  when  the  population  density  Is  high  and  when 
the  plants  are  healthy  and  have  not  been  attacked  by  Insects  or 

(see  Chapter  3).  When  entangled  petioles  slide  against  each  other 
during  relative  notion,  friction  resists  chat  notion.  As  a water 

liquid  exists  between  these  petioles,  because  leaves  are  vet.  If  the 
resulting  fluid  flow  Is  lanlnar.  the  following  di 
sliding  force  F : 


= B|_V„  (5.5) 


hyacinth  oat  with  respect  to  a spatial  reference,  and  B.  Is  the 
interactive  plant  friction  coefficient,  which  ia  proportional  to  the 

thickness  of  the  tiln  (Close  and  Fredrick,  1978),  Faraneter  B|^  was 
evaluated  through  slnulation, 

set  of  output  equations  is  derived  fron  the  state  variables. 

Mode  equations.  The  system  model  has  two  nodes:  therefore,  two 
vertex  equations  describe  force  interactions  at  these  nodes.  Applying 

zero.  The  node  equations  are 


2 -fc  * '"k  ' ° fS.7) 

State  variable  eeuatlona.  A set  of  state  variables  describes 
ccopletely  the  behavior  of  a system  In  response  to  its  Inputs. 

values  of  the  inputs  for  all  t 2 to  is  sufficient  for  evaluation  of 

Federlck,  1978).  One  state  varlble  is  usually  associated  with  each 
energy  storage  element  (e.g.  Inertia  and  spring  elements).  For  this 
study,  the  state  variables  are  Xg^.  V,  and  V^.  The  former 
pertains  to  cantilever  deflection,  and  variables  and 

correspond  respectivity  to  horizontal  velocities  of  the  front  and 


To  derive  state  variable  equations, 
of  eacli  state  variable  is  fomulated  in  tents  of  other  state 
variables,  tine,  and  Inputs  by  using  nodal  equations  and  component 
descriptions.  The  description  used  for  the  rate  of  change  of  Xg^ 
derived  below,  and  it  is  in  terns  of  input  variable  Vg  and  stati 
varloble  V,: 


16-8) 


The  state  variable  equation  for  V,  was  derived  in  the  followli% 

component  description  of  the  Inertia  of  the  enclosure  device  (Equation 

8.0): 


nodal  equation  5.8,  and  substituted  into  Equation  5.9: 


Finally,  the  second  state  variable  equation  was  obtained  after 


>.11) 


variable  aquation  {Equation  5.! 


^ ' V " ■ " 

' [\  ( V.  - V,  )]  (5.12) 

Output  eouationa.  Hat  conpaatlon.  tovlng  force,  sliding  contac 

and  they  are  calculated  at  any  tine,  t,  fron  the  state  variables. 
Slldltg  contact  and  drag  force  are  conputed  fron  Equations  5.4  and 
5.5.  Hat  compaction  for  a mt  of  original  length  L is  conputed 


It  eonpaction  - L-  (Xq  -X,)  (5.13) 

r the  fol loving  equations  ai 


Hethod  of  Solution  aj 
Equations  5.8.  5.11  and  5.11 

water  hyacinth  towing  system.  The  system  was  t 
to  reduce  the  towing  force,  study  mechanisms  oJ 


if  which  Is  nonlinear. 


ey  represent 


equaled  L.  Further,  at  t 2 to.  velocity  Vq  wae  a function  of  time, 
and  aystea  variables  and  outputs  were  evaluated  numerically  until  the 
system  reached  steady  state  or  until  the  derivatives  of  the  state 

The  fourth-order  Runge-Kutta  fonsula  was  used  to  solve  the 
equations  In  a "problem  solver"  program  written  in  Turbo-l^scal 
{Appendix  A)  The  program  utilizes  a published  version  of  Runge-Kutta 

seconds  for  an  Integration  interval  of  15  seconds. 

To  estimate  the  parameters  and  C^.  the  following  scheme  was 

adopted.  The  program  accepted  eatinmtes  for  parameters  and  C^, 

solved  the  differential  equations,  and  returned  values  of  towing  force 

parameter  estimation  was  repeated  until  correlation  between  calculated 
and  experimental  data  was  better  than  SOX,  Towing  force  on  oats  like 
those  described  In  towing  experiments  d through  7 {Table  d,0)  was 
calculated.  Plants  in  the  sacs  ranged  from  very  small  to  big. 


s asfiutDDd  CO  equal  the  experiinenclat  value 
It  do  not  vary  appreciably  In  a water  hyacinth  neit.  Friction 

hyacinrbs,  filn  thlclcnesa  and  viscosity  do  not  vary,  and  contact  area 

Reaults  of  i»«r««-t.r  Estlsatlon  and  Model  Anavlsls. 

Correlation  coefficience  between  nodal  and 

exporlnental  and  calculated  towing  force  ccnpared  at  Cels.  This 

Expressions  describing  B,  and  C.  selected  on  the  basis  of 
these  correlations  {Table  5.0)  were  found  to  depend  on  plant  size  and 

was  found  to  vary  with  Xj  until  steady  conditions  were  roaebed.  and 


{5.16) 


s necessary  b 


iss  sUbmii  It  would  no 

that  was  comprised  of  large  plants  and  large 
the  drag  eoeflclant  just  as  a mat  starts 


roots  and  large  rhizomes  would  effect  Ic.  In  addition,  D was  not 
affected  by  aat  geometry  or  cowing  velocity. 


Table  &.1  Values  of 


The  values  of  determined  for  each  mac  modelled  i 

in  Table  5.1.  It  was  affected  by  Che  type  ai 
that  comprised  a mac,  but  it  was  unaffected  1 
speed.  The  sac  with  the  lowest  value  of 
very  snail  plants.  The  friction  factor  of  rats  of  larger  plants 
varied  frmn  IdO  Co  ITS.  This  range  of  values  could  be  due  to 
connectivity,  Ssmll  plants  were  connecced  primarily  by  stolons,  and 
larger  plants  were  connected  primarily  by  overlapping  leaves. 

Effect  of  cowlne  velocity.  Macs  were  accelerated  ec^rlmentally 
in  one  of  two  ways.  For  cowing  velocltlos  less  than  0.25  m/s,  the 
acceleration  resembled  a step  velocity  input,  and  for  cowing 
velocities  ^ 0.2S  n/s,  the  macs  wore  accelerated  Incrementally  to 
offset  possible  detrimetitsl  inertial  effects  on  Che  testing  equipment. 


single  step  velocity 


Input.  When  a imt  ms  towed,  towing  force  Incrsesed  very  quickly  and 

Independently  of  doc  area  at  apprcxlsBtely  0.1  s,  because  Inertia  due 
to  the  mess  of  the  enclosure  Is  the  principal  active  cooponent  on  the 

frequency  of  the  spring  - enclosure  systea  was  5 hz.  As  speed 
increased  and  mat  length  Increased,  it  took  longer  for  the  force  to 
deeroase  to  a steady  value  (Figures  S.l  and  5.2).  For  a sanple  size 
2.4d  a long  by  1.22  to  wide,  tine  to  reach  steady  conditions  ranged 
fron  2 s for  a velocity  of  0.9  n/a  to  10  a for  a velocity  of  0.2  ntfs. 

Sliding  friction  (Equation  5.d  ] and  viscous  friction  (Equation 
5.5)  are  responsible  for  danqilng  (Figure  5-3).  Viscous  friction 

It  quickly  becomes  the  principal  active  component. 

Mat  compaction  as  measured  at  steady  state,  increased  linearly 
with  tawing  velocity.  When  B.  equaled  30,  compaction  of  a 2.dd  m 
long  X 1.22  m wide  m iDt  increased  with  towing  velocity  according  to 
this  expression: 


To  offset  the  effects  of  inertia  of  the  enclosure,  mats  were 
accelerated  incrementally  whenover  the  desired  final  towing  velocity 
was  greater  than  0.25  m/s.  Increments  varied  according  to  the 
discretion  of  the  trolley  operator.  Velocity  inputs  to  the  problem 
solver  were  changed  Incrementally  with  tine  in  an  attempt  to  mimic 
these  experiments.  Two  simulations  were  done,  one  involved  a eat  with 
very  small  plants;  the  other  involved  isedlun-sized  plants  typical  of  a 
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«ater  creatinenc  system.  The  former  is  presented  In  Figure  S.4. 
the.t  every  tine  velocity  was  incremented,  there  would  bo  a peak  In  the 


not  great  enough  to 


to  the  sparseness  of  ekperlmental  data  in  th' 

On  the  other  hand,  a mat  comprised  of  larger  plants  had  friction 

experimental  data:  although,  the  model  accurately  predicted  the 
general  Increase  in  towing  force.  This  discrepancy  problably  was 
related  to  the  practical  impossibility  of  a step  velocity  input. 

The  advantage  of  reducing  Inertial  effects  through  a slow 
increase  in  towing  velocity  is  that  the  eaxinun  towing  force  can  be 
reduced  by  a factor  of  four  (Figures  5.S  and  S.6).  The  "problem 

other  kinds  of  velocity  inputs.  Different  ramp  velocity  inputs  were 
modelled  (Figure  5.7),  and  results  indicate  that  inertial  forces  could 


negligible  for  snail  accelerations. 
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tovlng  force  is  evident  in  Figure  5.1 
density,  20S/a)  have  small  rhizome. 

cooqsared  to  the  expected  total  area 


lie  effect  plant  size  has  on 

and,  although  the  total  numi 
projected  area  in  flow  is  sf 


Therefore,  as  depicted  In  Figure  5,8,  the  viscous  force  at  low  towing 
velocities  on  small  rhizomes  was  not  large  enough  to  dampen  quickly 
the  oscilatory  response  of  the  ccsnblned  eprlng-mss  components.  In 

Also,  Initial  oscilatory  response  of  spring-mass  components  takes 

mat  decreased,  sliding  force  contributed  a greater  proportion  of  the 

5.10  and  5.11.  Table  5.2).  Compaction  of  a 1.22  m long  by  0.51  re  wide 
mat  was  approxleretely  half  the  compaction  of  a 1.22  re  long  by  1.22  m 
wide  sat.  On  the  other  hand,  as  the  Initial  length  of  the  mat 
decreased,  net  compaction  (X)  Increased  slightly. 


•d  plants  averaging  61.0 
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the  body  tips.  While  A water  hyacinth  Is  vertical  due  to  support  by 
other  plants  or  by  floats,  Its  centroid  of  displaced  volume  Is 


lies  above  the  cshtrlod  of  displaced  voluise.  These  ci 
water  hyacinths  unstable  vertically  (Gerhart  and  Gross.  1985).  Pla] 

If  the  plants  are  tipped  so  that  I 

buoyant  force  and  centerline  of  the  volume  of  the  object. 

The  buoyant  force  of  a eat  was  described  in  chapter  3.  Buoyant 
force  per  unit  sat  area  depends  on  areal  bionas  density  and  on  tbe 
physical  condition  of  the  plants.  If  plants  are  well  connected,  force 

the  support  structure  of  a mac.  Hats  cowed  In  this  study  encored  Che 
under  occurred  ranged  between  0.4  and  0.6  is/s.  Calculations 


submerge  them  (e.g.  a 3 or  imot  required  over  500  N to  submerge  it), 
plants  and  connectivities  become  unstable  at  different  values  of 

of  larger  plants  become  unstable  when  compaction  was  approximately 


One  of  the  assumptions  made  1] 


e Independent  of  sample  size  if  the  mat 
saiac  size  and  areal  density. 

consistently  "roll  under"  at  lower  compaction,  the  model  was  changed 
are  In  contact  within  the  mat  Co  prevent  further  compaction).  To 

nature  of  rhizomes  In  contact.  It  is  described  as  follows: 


Ky  = 10  xlO^^  Pa  - The  elastic  constant  Is  the  elastic  sw>du 

= [K  Xo.-lOOO  Hp  V,  IV.  |-  Bj^(V,-V,) 

- ( L e C X.  - X,  )-  Dp/M^  (6.17) 


Resolts  showed  that  when  compaction  equaled  D . towing  forces  per 


sufficient  CO  subcoerge 


Decs,  prcpvldod  Chat  the  plants  contact  at  an  angle,  and  a vertical 


.e  dlnensiana]  problens. 


-Hats  of  water  hyacinths  towed  within  a recconglar  enclosure  were 
friction  eletwent,  an  elastic  element  For  the  enclosure,  and  an  elastic 


different  sizes  and  types  of  nets  and  velocity  Inputs. 
- A DCthod  to  predict  when  a nat  becomes  unstable 


‘on  Che  opposlce  end.  Cblculaced  compressive  force  was 
I evpcrlnencal  data.  This  inforiDaclon  can  be  used  Co 


Conceolual  Hodel. 

The  conceptual  model  of  a vacer  hyacinth  mat  of  length  L and 

represent  Che  water  hyacinth  mac.  The  compressive  frame  represents 
Che  test  equlimaent  used  to  compact  mats,  so  data  obtained  from 
experiments  could  be  used  to  validate  Che  sat  model.  The  sac  was 


situated  between  too  bars,  a front  pusher  bar  and  a fixed  rear  bar. 
The  front  bar  pushing  the  eat  was  considered  rigid;  therefore,  only 


Initially  the  system  is  at  rest.  The  distance  from  node  0 to 


IIS 


velocity  Vg,  it  coB^resses  the  net  a^inst  the  back  bar,  and  the 
shortens.  Compaction  continues  until  the  plants  "roll  under." 

The  element  descriptions  were  given  In  Chapter  5.  Newton's 
law  applied  at  node  1 produces  Che  following  equation: 


= K,  ( X.  - L ) 


c force  of  plants  against 


bar.  and  it  equals  zero  when  X,  - L,  Nest  length  of  the  m: 

X lO'  N/n,  which  equals  approxineteiy  the  Inverse  of 
modulus  of  elasticity  of  Che  two  contacting  bodies  (rhizome  and  steel 
rear  bar)  multiplied  by  eat  width  V.  The  state  variable  equations 


, ) - ( X,  - L ) ]/K_,  (6.5) 


al  conrpreB9lve  force.  Compressive  force  was  calculated  at  node 
«i  input  variables  and  computed  state  variables  for  time  t ^ tg 

‘ "a  H"  * "p  *r  *^d  * ®L 


t for  a ramp  velocity  input. 


Hat  oonpection  *as  calculated  on  the  basis  of  the  change  in 
Initial  areal  standing  densities  of  mats.  Specific  area  A is  defined 


*,  = {1.  - Xo)  » 

where  Xg  was  obtained  after  the  integration  of  Vg  - . Q. 

Method  of  Solution  pH  ParHITter  Estimation. 

Equations  €.4  and  6.5  represent  a system  of  coupled  first-order 
differential  equations.  Together  with  the  output  equations,  they 
represent  ecsopaction  of  a water  hyacinth  nat.  The  equations  were 
solved  to  determine  the  behavior  of  a mat  subjected  to  different 
compaction  velocities,  mechanisms  of  "roll  under",  and  parameters 


end  Cj  of  unaceady  flow  oondltlona.  Tba  latter  were  studied  b< 
To  faclllate  the  estiiBtlon  process  Shd  to  solve  the  syatei 

conditions  were  V,  = 0,  XisL  and  Xg  s 0.  Input  velocity  wa 
function  of  tlue,  and  It  was  specified  to  be  a ramp  In  the  firs 


were  chosen  because  they  modelled  accurately  the  experimental 
was  used  to  attempt  to  produce  a better  correlation  between 


experimental  and  calculated  eoaipresslve  force. 


Initially  varied  according  to  a 


(6.8) 


EicBtrlJOTIal  Conor 


modified  automotive  winch.  With  these  devices  four  different. 


through  a raJte  support  carriage  equipped  with  load  cells  to  measure 


position  sensing  potentiometer  were  recorded  on  a Campbell  Scientific 
an  data  logger  equipped  with  "Burst  mode  ' software  which  allowed  a 


sampling  Interval  of  0.2  s. 

Tests  wore  conducted  on  plants  at  the  following  locations:  IFAS 


{Figure  6.2}.  seven  plants  were  characterized  by  shoot  height,  root 


plants 


weighed,  and  the  weight  was  divided  by 


Hat  buoyanc 


CoBpression  frame  was  equipped  with  a eliding  bar  and  load 
cells  located  at  both  coda  of  the  bar.  As  the  bar  moved 
toward  the  flrod  end  bar,  the  not  compacted, 


Figure  6.2  A water  hyacinth  mat  was  being  conpacced. 


I 15  dlfferenc  in 


of  conpacclon  velocity,  coimeotivlcy.  era 


(Nervations 


well  connected 


0.22  Elckly 


Before  a compaction  emperlnent  started,  a grid  of  nylon  cord  iras 

under  ~ occurred  because  It  begins  vhen  grid  lines  bunch  closely 
together.  All  experiments  mere  videotaped.  After  the  tapes  mere 
viewed,  Che  position  of  the  moving  ralca  at  "roll  under"  mas  recorded. 
Results  and  Discussion. 

force  compaction  behavior,  "roll  under"  and  klnenatlcs.  (lompaetion 
experiraonts  Indicated  that  ae  compaction  velocity  increased,  the 

velocity  {<  0,15  m/s],  the  relation  between  coo^ressive  force  and 
specific  area  mas  approxin&tely  linear  up  to  "roll  under".  Ac  higher 
velocity  (>  0.25  m/s],  the  relation  was  obviously  nonlinear.  The 


Paraieter  deterilnallun.  Paranelars  and  Bj^  varied  as  a 


(Close  and  Frederick.  1978).  Since  water  hyacinth  contact  area 


Table  6.1  Contact  friction  coefficient  and  ocher  variables 
used  In  the  Blnjlacione  of  nac  compaction. 


Exp. 


Compaction 

velocity 


Friction  coefficient  8^^  decreased  with  compaction  velocity  and 
connectivity  (Table  6.1).  Plants  in  experiments  3 and  S »cre  of  low 

highly  connected.  Friction  coefflclancs  of  different  mats  compacted 
at  Che  some  velocity  was  always  significantly  tower  In  uses  of  lower 
connectivity.  Friction  eoefflelent  of  well  connected  oats  decreased 
linearly  with  cooipacclon  velocity,  and  Che  regression  relation  follows 
(r®  = 0.96) 


ta  points  were  available 
'W  connectivl  ty;  however, 
>se  of  highly  connected  IS! 


velocities  exceeding  D.26  m/: 
An  explanation  for  t 


a Function  of  velocity  can  be  found, 
d linearly  with  towing  velocity 
s Inpllea  that  a imt  reaches  a maxlmun 

"roll  under"  in  a nat  Is  acbelved  at  loser  compaction  velocity  by 
restralnclng  It  and  forcing  it  to  compoct.  Therefore,  extra  friction 
force  la  required  at  low  cooipactlon  velocities  to  ccsoixici  a oci  until 

compaction  velocity,  the  fraction  of  total  compressive  force  due  to 
contact  friction  was  great:  at  higher  velocities.  It  was  snail. 

Graphs  of  experimental  data  showed  that  at  low  compaction  velocity  the 

(Figures  6-5  - 5.9). 

compressive  force  in  a nat  that  consisted  initially  af  amny  open 
spsacea  (Figure  6.9).  Connectivity  Is  considered  to  be  zero  in  these 
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6.6},  because  the  mess  element  representing  Che  mass  of  Che 

Ccmpecclon  (X)  at  ~roll  under"  varied  according  to  plant 
population  and  plant  length.  Generally,  It  decreased  as  population 

betweeo  20  and  40X  conpactlon.  which  was  higher  than  the  values  In 


dirricuic  fo 


. "Roll  under"  occurred  at  the  front  end  of  a mat  during 
but  Inside  the  sat  during  compaction  tests.  Therefore, 
naturally  supportive  in  coegwiccion,  making  it  more 

was  large  enough  Co  sufaewrge  the  mat,  provided  Che  plants 

large  plants,  "roll  under"  was  usually  accompanied  by  a 


-A  matheimtlcal  model  was  developed  to  predict  Che  response  of  a 
water  hyaclath  nat  to  forced  compaction.  Initially  unknown 
partUDecers,  and  3^,  were  estimated  through  comparison  of 

calculated  and  eKporlmental  compressive  force..  Friction  coeffletenc 
varied  with  compaction  velocity,  compaction,  and  connectivity.  Drag 
coefficient  varied  with  Che  change  in  length  of  a obc. 

-The  model  for  mac  compaction  accurately  predicted  compressive 
force.  "Roll  under"  con  be  predicted  on  a basis  of  corapaoclon,  or 
buoyancy . 

-Cdnpressive  force  wae  dominated  by  drag  at  velocity  exceeding 
0.26  m/s  and  ^ contact  force  at  velocity  below  0.16  n/s. 


CHAPTER  VII 

CntCLUSIC»S  AND  REOWODATICWS 
FOR  FORTHEER  RESEARCH 


40X  (Coapuced  (lengch/1-}100) , 
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obtaiiMd.  Poisson’s  ratio  and  Young's  siodulus  sore  dccemlned, 
respoccivlty,  to  be  0.47  and  9.9  MPa. 

hyacinth  interactions  for  two  different  harvesting  operations: 
rectar^lar  itat  towing  and  forced  coepaetlon.  Hass  and  elastic 
elenents  of  test  oBchlnery.  petiole  contact  friction  elenent,  moss 

Their  descriptions  and  Interconnections  were  based  on  experinontal 
observations  and  physical  and  hydronwtchonical  properitles.  Tho 
equations  derived  fron  the  conceptual  aodel  were  solved  nusierlcally 
using  a Rur^e-Kutta  integration  method  to  study  system  behavior  in 
response  to  velocity  input.  Output  variables  such  as  towing  or 
compaction  force  and  sett  cct^ction  were  predicted  over  time  with  good 
ith  cjcperliiiontol  data.  Hat  Instabtlty  or  "roll  under" 

1 for  the  design  of  be 


harvesting  schemes  and. 

-Effect  of  different  acceleration  patterns. 

-Power  and  energy  requirements  of  different  designs. 

The  design  of  equipment,  such  as  conveyors,  rakes,  removal 
devices,  and  restralntlng  barriers,  could  be  based  on  the  following 

-Slowing  accelerating  a mt  and  moving  machine  parts  to  reach  a 
constant  velocity  saves  energy.  Because  tho  highest  foroes  in 


Che  syscens  studied  were  due  CO  Inertia,  plants  should  be  kept 
moving  at  constant  velocity. 

-Drag  force  is  reduced  by  cowing  long,  narrow  eats. 

-"Roll  under"  is  prevented  by  maintaining  approplate  plant 
density  and  operating  speed. 

Additional  research  should  be  done  for  large  scale  verification 
and  application  of  Che  model.  Moreover,  model  parameters  connectivity 
and  petiole  friction  coefficient  could  be  Investigated  to  determine  on 
empirical  relation  between  the  two  parameters.  Dependency  of  this 
coefficient  on  other  plant  characteristics  could  be  clarified,  too. 
Also,  because  C.  was  determined  only  for  recCanglar  mats,  drag 

program  with  input  and  output  screens  Is  needed  for  user  convenience. 

Finally,  the  movement  only  of  a mat  and  moving  machine  parts. was 
investigated.  Study  of  water  flow  under  a moving  mat,  which  Is  a very 
complicacod  problem,  could  describe  mat  movement  due  to  negative 
acceleration. 
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> v[2];5;3.’  •.  dl8Witfe-6;3.  ' v[3];6;3,'  v[4]:S: 

> -.O.O.'  •.{length-v[3]«[0):6:3 

> .•  •,vel:5:3.'  • , trlal:6:2) : 

> h :=  (x2-xl)/nstep: 

> derlvsfx.v, delta, dv); 

> rk4(v.dv.nvar,x,h. delta, vout); 

> wricelnC'peuse  in  routine  R10UNB‘): 

> ericeln( 'stepsize  to  snail'};  readln 


if  ( x>e2,2S)  and  (x  < 3.S)  then  begin 


if  { x>=3,5)  and  (x  < 5.50)  then  begin 


> and  (v[4]  >v[2])  or  ( constant  Mength) 

> then  distancel  1 ;slengch  else 

> dlstancell;ev[l]+conscant; 

> trial  :=Hp«Arx0.4245xexp(0.933*nridth/diatQiiceU)aahs(v[2]}«v[2] 

> (output  variables,  remove  the  brackets  around  fout  to  write  to 


>r  Compaction  Oiora< 


M<m.vu..i 


R.  L. 


K^16-? 


iS-~= 
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